Experimental study of rj c decays into vector-vector final states 
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Using 58 million J/ip events accumulated with the BES II detector, branching fractions for the 
rj c — > VV decays, rj c — » pp, rj c — > K* K , r) c — * uju), and rj c — > (f><f>, are measured. The branching 
fractions are Br(-q c -> pp) = (1.25±0.37±0.51) x 1(T 2 , Br^ -> K*K*) = (10.4± 2.6±4.3) x 1(T 3 , 
Br(ri c -» ujuj) < 6.3 x 10~ 3 (90 % CL) , and Br(r) c -» H) = (2.5 ± 0.5 ± 0.9) x 10~ 3 . The process 
rj c — > LO(f) is also searched for, and the 90% CL upper limit Br(r/ C — > uj4>) < 1-3 x 10~ 3 is obtained 
for this process. 
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PACS numbers: 13.25.Gv, 14.40.Gx, 13.40.Hq 

I. INTRODUCTION 

It has been known for many years that r\ c decays into vector-vector (VV) mesons, even though this process 
is forbidden by hadron helicity conservation (HHC) assuming collinear valence quark configurations dominate. Three 
different Bethe-Salpeter wave functions (Gauss, Exponent, Power) have been used to calculate Br(ry c — > VV) [2]; the 
calculated results are lower than the experimental ones from the PDG [1] by two or three orders of magnitude, as 
shown in Table I. New theoretical predictions are needed. 

TABLE I: Theoretical predictions for Br(rj c — > VV) compared with PDG results [1]. Gauss, Exponent, and Power refer to 
three different Bethe-Salpeter wave functions used to calculate Br(j) c — ► VV) [2]. 
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Exponent 
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Br(j] c -> pp) 
Br{r) c -> K*K*) 
Br(j] c -> 4><f>) 


2.3 x 10" b 
2.8 x 1(T 5 
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8.7 x KT 5 
8.6 x 1(T 5 
1.6 x 1CT 5 


2.8 x 10" 4 
2.8 x 1(T 4 
5.0 x 1CT 5 


(2.6 ±0.9) x 10^ 
(8.5 ±3.1) x 10 -3 
(2.6 ±0.9) x 10" 3 



Branching fractions and reduced branching fractions, BRr, compared with SU(3) predictions [4] are shown in 
Table II. Here the reduced branching fraction is defined as BRr = BR/(P V ), where Py is the momentum of the vector 
meson in the r\ c center-of-mass system. The experimental results are not inconsistent with SU(3) expectations, but 
due to the large errors no meaningful investigation of possible SU(3) symmetry breaking patterns can be performed. 

TABLE II: Comparison of PDG Br(r/ C — > VV) results with predictions of SU(3) symmetry. Ratio of BRr is defined as 
BRft(rj c — ► VV) I ' BRr(t) c — ► (j)(j>), and this can be compared with the SU(3) expectation, SU(3). 





BR(10"- 5 ) 


br r (w-->) 


Ratio of BRr 


SU(3) 


Br{r) c -> pp) 


26 ±9 


12.6 ±4.4 


6.3 ±2.3 


3 


Br( Vc ~^K*K*) 


8.5 ±3.1 


5± 1.7 


2.5 ±0.9 
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2.6 ±0.9 


2 ±0.7 


1±0.4 
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Exclusive decays of r\ c into vector meson pairs have also been investigated in the framework of the 3 Pq quark- 
creation model [3]. The width of r/ c — » p a p° in PDG is significantly larger than the model prediction. The branching 
fraction of rj c — > uxj) is very important to determine the mechanism of r\ c — > VV decay [5]. In this paper, J/ip — > jr) c 
decays from a sample of 58 million J/ip decays obtained with the BESII detector are used to determine Br(j] c —* pp), 
Br(j] c —> K*K ), Br{r\ c — > lulu), Br{r\ c — > 4>(f>), and Br{r\ c — > LU(f>). 

II. BES DETECTOR 

BES is a conventional solenoidal magnet detector [6, 7]. A 12-layer vertex chamber (VTC) surrounding the 
beam pipe provides trigger and trajectory information. A forty-layer main drift chamber (MDC), located radially 
outside the VTC, provides trajectory and energy loss (dE/dx) information for charged tracks over 85% of the total 
solid angle with a momentum resolution of <J p /p — 0.0178^/1 + p 2 (p in GeV/c) and a dE/dx resolution for hadron 
tracks of <~ 8%. An array of 48 scintillation counters surrounding the MDC measures the time-of-flight (TOF) of 
charged tracks with a resolution of ~ 200 ps for hadrons. Radially outside the TOF system is a 12 radiation length, 
lead-gas barrel shower counter (BSC). This measures the energies of electrons and photons over ~ 80% of the total 
solid angle with an energy resolution of <te/E — 2\%j\[E (E in GeV). Outside the solenoidal coil, which provides a 
0.4 Tesla magnetic field over the tracking volume, is an iron flux return that is instrumented with three double layers 
of proportional counters that identify muons of momentum greater than 0.5 GeV/c. 

A Gcant3 based Monte Carlo, SIMBES [8], which simulates the detector response, including interactions of 
secondary particles in the detector material, is used in this analysis. Reasonable agreement between data and Monte 
Carlo simulation is observed in various channels tested, including e+e~ — > (j)e + e~, e + e~ — > {"f)pp, J/tp — > pp, 
J/tp — > pir, and ip(2S) — > Tv + n~J/ip, Jjip — > For Jjip — > jr] c , decays are generated with an angular distribution 

of 1 ± cos 2 9, where is the angle between the e + and rj c in the laboratory. 
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III. Vc — p°p° 
A. Event selection 



First a J/ip — > 77r+7r~7r + 7r~ sample is selected. Events are required to have four good charged tracks and 
one or more photon candidates. A good track, reconstructed from hits in the MDC, must be well fitted to a helix 
originating from the interaction point; have a polar angle, 9, with | cos#| < 0.8; and a transverse momentum greater 
than 60 McV/c. TOF and dE/dx information are combined to form particle identification probabilities for the pion, 
kaon, and proton hypotheses. At least three tracks must be identified as pions. To reduce the number of spurious low 
energy photons produced by secondary hadronic interactions, photon candidates must have a minimum energy of 50 
MeV and be outside a cone with a half-angle of 15° around any charged track. 

To get improved momentum resolution and to remove backgrounds, events are kinematically fitted to the 
J/ip — > 77r + 7r~7T + 7r~ hypothesis, using all photon candidates. The fit with the highest probability is selected, and the 
X 2 of the four constraint (4C) fit is required to be less than 20. 

For J/ip — > 77r + 7r~7r + 7r~ decay, a major source of background is from J/ip — > 7r 7r + 7r~7r + 7r~. To remove events 

containing a 7r°, when there are multiple photons, |to 7 i 7 2 — m^ol > 60 MeV/c 2 is required if P m iss is near the plane 
of the two photons, 7I and 72, i.e. \P m iss ■ (ffi x ^72)! < 0.15. Here P m is.s is the unit vector of the missing momentum 
of all charged tracks; f 7 i and f 7 2 are unit vectors in the 7I and 72 directions, respectively; and to 7 i 7 2 is the invariant 
mass of 7I and 72. Two additional requirements, x 2 (J/V' - * 77r + 7r~7r + 7r~) < x 2 {J/tp — > 777r + 7r~7r + 7r~) and P i7 < 
0.0015 (GeV/c) 2 , are used to further remove J/ip — > 7r 7r + 7r~7r + 7r~ background. P tj is the transverse momentum of 
the 7r + 7r~7r + 7r~ system with respect to the photon. Finally, the requirement \U m iss \ = \E m iss — cP m iss \ < 0.07 GeV/c 
is used to reject events with multiple photons or charged kaons; here, E m i SS and P m iss are, respectively, the missing 
energy and missing momentum calculated using only the charged particles, which are assumed to be pions. 

There are other possible backgrounds such as J/ip — > W7r + 7r~, J/ip — > jK®K®, and J/ip — > jK®K ± n T . 
The J/ip — > aj7r + 7r~ background is suppressed by the requirement \M^+ n - n o — M u \ > 40 MeV/c 2 , where the tt° in 
7r°7r + 7T~7r + 7r~ is associated to the missing momentum and energy, determined using only the charged tracks. To 
remove the J/ip — > jK®K® background, \M n + w - — M K o | > 25 MeV/c 2 is required for both ir + ir~ pairs. To remove 
the background from J/ip — > -fK®K ± iT T , events are rejected if x 2 (J/?/> — > jK ± ir T n + ir~) < x 2 (J/ip — » 77r + 7r~7r + 7r~) 
when \Mk+k- - M K « \ < 25 MeV/c 2 . 

Figure 1 shows the scatter plot of M^ttZ versus M vr i vr4 for surviving events with M 7r + 7r - 7r + 7r - > 2.8 GeV/c 2 , 
where nl (tt3) and 7r2 (7r4) are the 7r + and 7r~ with the higher (lower) momentum. There are clear signals near (M„o, 
M p o) and {Mf 2 n270)i -^/ 2 (i270))- The decay J/ip — > 7/2/2, including ?7 C — > /2/2, has been studied previously using 
the J/V> events collected by BES II [9]. 
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FIG. 1: Scatterplot of M W 2n3 versus M n i n 4. Also shown are the signal and sideband background boxes used in this analysis. 



The number of r\ c —> p°p° events and the corresponding backgrounds are estimated from the boxes in the 
scatter plot of tt + tt~ versus 7r + 7r~ invariant masses, as shown in Fig. 1. The signal region is shown as the square box 
22 at (0.77, 0.77) GeV/c 2 with a width of 240 McV/c 2 . The main backgrounds are r\ c — > pir+ir^ and r\ c — > 7r+7r~7r+7r~. 
The T]c — * pn + ir~ background produces the horizontal and vertical bands at m p , and the i] c — » 7r + 7r _ 7r + 7r _ background 
produces the uniform background in the m^+^- versus m^+ v - scatter plot of Fig. 1. These backgrounds are estimated 
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FIG. 2: Fitting r\ c in the M n +^-^+ v - spectrum of events in the signal box 22 in Fig. 1. The fitted number of r\ c events is 
N 3i9 = 106.9 ± 18.2. 



using the horizontal and vertical sideband background boxes, which are taken 50 MeV/c 2 away from the signal box 
and denoted as 12, 21, 32, and 23, and the two diagonal boxes, denoted as 31 and 13, in Fig. 1. 




FIG. 3: Fitting r\ c in the M n + w - w + n - spectrum of four horizontal and vertical sideband boxes (12,21,32,23). The fitted number 
of rj c events is N sidl = 83.3 ± 21.5. 



B. Results 



Shown in Figs. 2, 3, and 4 are the M 7r + 7r - 7r + 7r - fitting results for p°p° signal and sideband regions. Binned 
maximum likelihood fits are performed with a Breit-Wigner convoluted with a Gaussian resolution function for the 
rjc and a first order polynomial for background. The number of r/ c — > p°p° events is obtained with the formula 
N obs = N st 9 _ N sidij 2 + N sld2 /2. Here N sia , N sU1 , and N sld2 are the fitted number of r] c events in the Af 7r+7r - 7r+7r - 
signal box and sideband boxes(See Figs. 2, 3, 4). So N obs = (106.9±18.2)-(83.3±21.5)/2+(15.8±8.9)/2 = 73.1±21.6. 

The r] c — > p°p° branching fraction is determined from 

Br{llc p ° p0) = iVW^) 

= (4.15 ± 1.23) x 10" 3 , 
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FIG. 4: Fitting r\ c in the M 7T + 7T - 7T + 7I - spectrum of two diagonal boxes (31,13). The fitted number of rj c events is N' 
15.8 ±8.9. 



where e — 2.35% is the efficiency, corrected for sideband subtraction, determined from Monte Carlo simulation. 
Nj/.,p = 57.7 x 10 6 is the total number of J/ip events collected by BESII [10], and Br(J/ip -> jr] c ) = 0.013 is the 
branching fraction for J/ip — > ^r\ c [1]. Correcting for charged decay modes, Br{r\ c — > pp) = (1.25 ± 0.37) x 10~ 2 is 
obtained. 



IV. J] c -f K* K* 

J/ip — > 7^ c — > jK*°K*° — » 77r + 7r~_ft' + _ft'~ is studied to determine the r/ c — > K*K* branching fraction. 

A. Event Selection 

The selection criteria for charged tracks and photons are the same as for the r\ c — > analysis. At least one 
good photon is required. TOF and dE/dx information are combined to form particle identification probabilities for 
the pion, kaon, and proton hypotheses for each track, and at least two tracks must be identified as kaons. The same 
method as used in the r\ c — > p°p° analysis is used to remove n° background. 

Some additional requirements are applied. Combined probabilities are formed from the particle identification 
probabilities of all tracks and the probabilities determined from the 4C kinematic fit to the J/ip — ► 77r + 7r~ K + K~ 
hypothesis for all possible combinations. The fit with the highest probability is selected, and the x 2 of the 4C 
kinematic fit is required to be less than 10. 

To further remove backgrounds containing 7r°s, two requirements, X 2 {J ^ 77r + 7r~ K + K~) < 
X 2 (J/'0 — ► 777r+7T~ K + K~) and P% < 0.003( GeV/c) 2 , are used. Pt 7 is the transverse momentum of the tt + tv~ K + K~ 
system with respect to the photon. Furthermore, the requirement |f/ TO j ss | = \E miss — cP m i SS \ < 0.1 GeV/c is used to 
reject events with multiple photons; here, E m i SS and P m iss are, respectively, the missing energy and missing momen- 
tum calculated using only the charged particles, which are assumed to be tt + tt~ K + K~ . x 2 {J ' li> ir + ir~ K + K~) > 20 
is required to remove J/ip — ► tt + tt~ K + K~ background. 

There are other possible backgrounds such as J/ip — ► uK + K~, J/ip — » jK °K ± ir T , etc. The J/ip — > luK + K~ 
background is suppressed with the requirement {M^+^-^o — M u \ > 40 MeV/c 2 , where the ir° in ir°ir + ir~ K + K~ is 
associated to the missing momentum and energy, determined using only the charged tracks. To remove the background 
from J/ip — ► jK°K ± n T , events are rejected if x 2 (J/'0 jK ± ir T n + Tr~) < \ 2 {J/ip — * jtt + tt~ K + K~) when \M v + n - — 
M K o\ < 25 MeV/c 2 . To remove backgrounds containing </>, \M K + K - — M$\ > 0.02 GeV/c 2 is required. 

A K*°K signal is seen in the scatter plot of M^+n- versus Mx- V +, shown in Fig. 5. The same method is 
used as for the r\ c — > p°p° analysis to estimate sideband background. The signal region is shown as a square box at 
(0.896,0.896) GeV/c 2 with a width of 160 MeV/c 2 . The main backgrounds are r] c -> K* Kn and i) c -» K+K-jt+tt-. 
The rj c — » K*Kn background appears as the horizontal and vertical bands at mjf., and the ?7 C — » A" + A" _ 7r + 7r _ 
background as the uniform background in the M K + W - versus M K -^+ scatter plot of Fig. 5. The backgrounds 
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are estimated from the sideband boxes, which are taken 40 MeV/c 2 away from the signal box and shown as four 
horizontal/ vertical boxes (denoted as 12,21,32,23) and four diagonal boxes (denoted as 31,13,11,33) in Fig. 5. The 
horizontal and vertical sideband boxes (12,21,32,23) allow the determination of the the backgrounds from the horizontal 
and vertical bands, while the diagonal boxes (13,31,11,33) allow the estimation of the smaller uniform background 
contribution. According to a Monte Carlo study, other -q c backgrounds such as rj c — > K + K~K + K~ and i] c — ► 
7r + 7r~7r + 7r~ do not survive the selection criteria. 
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FIG. 5: Scatter plot of M K+7r - versus M 



K-tv+ 



Signal and sideband boxes used in the analysis are shown in the scatter plot. 




2.8 2.85 2.9 2.95 3 

M 



KKjra 



FIG. 6: Fitting r\ c in the M n +„- K + K - spectrum of events from signal box 22. The fitted number of rj c events is N slg 
80.3 ±13.6. 



B. Results 

Shown in Figs. 6, 7, and 8 are the M n +„- K + K - fitting results for K*°K*° signal and sideband regions. 
Binned maximum likelihood fits are performed with a Breit-Wigner convoluted with a Gaussian resolution function 

for the rjc and a second order polynomial for background. The number of r\ c — > K*°K events are obtained from 
N obs = N sig _ N sidi/ 2 + N sld2 /A. Here N si a , N sidl and N sld2 are the fitted number of ^ events in the M^ +7! - K+K - 
signal and sideband boxes(See Figs. 6, 7, 8). So N obs = (80.3 ± 13.6) - (41.0 ± 10.2)/2 + (0.0 ± 0.1)/4 = 59.8 ± 14.5. 
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FIG. 7: Fitting r/ c in the M^+^-^+^f- spectrum of events from the four horizontal and vertical sideband boxes. The fitted 



number of r\ c events is N' 



idl 



41.0 ± 10.2. 




FIG. 8: Fitting ij c in the M n + n - K + K - spectrum of events from the four diagonal boxes (31,13,11,33). The fitted number of r\ c 
events is iV sid2 =0.0 ±0.1. 



The branching fraction of r) c — > K*°K* is given by 

ik* - A-»r») . _, u 

= (5.2 ± 1.3) x 10~ 3 , 

where e = 3.4% is the efficiency, corrected for sideband subtraction, determined from Monte Carlo simulation. 
Using Br(K*°(K* ) — > if^T) = 0.67 and correcting for charged decay modes, we obtain Br{r\ c — ► K*K*) = 
(10.4 ±2.6) x lO" 3 . 



V. ?7 C — ► 

J/^> — > 77? c — > 7^o) — > 77r°7r + 7r _ 7r°7r + 7r _ is studied to obtain the branching fraction of r] c — > uu. 
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A. Event Selection 



The charged track selection criteria and particle identification are the same as in the above two analyses. Here, 
four tracks must be identified as pions. The photon selection criteria are different for rj c — > lulu to get higher efficiency 
because there are five photons in this channel. To reduce the number of spurious low energy photons produced by 
secondary hadronic interactions, photon candidates must have a minimum energy of 50 MeV and be outside a cone 
with a half-angle of 6° around any charged track. At least five good photons are required. 

To get improved momentum resolution and to remove backgrounds, events are kinematically fitted to the 
J/ip — > 577r + 7r~7r + 7r~ (4C fit) and J/ip — > 77r 7r 7r + 7r~7r + 7r - (6C fit) hypotheses, for all combinations of photon 
candidates, xlci^/^ ~~ * 577r + 7r~7r + 7r~) < 20 and xicl-VV' ~> 77r 7r 7r + 7r~7r + 7r~) < 8 are required. To get cleaner 
J/ip -> 'jlulu candidates, \M^ - 0.135| < 0.035 GeV/c 2 and |M 2 - 0.135] < 0.035 GeV/c 2 are required. Here, the 
photons used for M\ and M 2 are determined from the best 6C fit to the J/ip — » 77r 7r°7r + 7r~7r + 7r~ hypothesis. 



B. Results 



There are four pairs of M^ +7r _ 7r0 versus M 2 +jr _ 7r0 combinations for every event. If one of the combinations 
satisfies \M* +n _ n0 -M u \ < 0.035 GeV/c 2 and \M% +n _ n0 -M u \ < 0.035 GeV/c 2 , the event is selected as a J/ip -» jlulu 
candidate. The M^o^+^-^o^+t^- distribution for J/ip — > ^lulu candidates is shown in Fig. 9. 




FIG. 9: The ij c fit of the M w o 7T + 7I - 7T a 7T + 1T - distribution of candidate events in the lulu region. 



Figure 9 shows the result of fitting rj c to the M v o v + v - v o v + v - distribution of candidate events in the lulu 
region. The shape of the rj c used in the fit is obtained from Monte Carlo simulation of rj c — > lulu. The fitted number of 
7] c events is 16.4 ± 9.4. The result has large errors and includes decays other than J/ip — > ^lulu decay. Therefore only 
the upper limit on Br(rj c — > lulu) is determined. At the 90% confidence level, the upper limit on N Vc is N upper = 23. 
The upper limit on Br(?7 c — ► lulu) at the 90% C.L. is determined from 

]\J upper J e 

Br{j] c — > lulu) < 



Nj/^ ■ Br(J/ip — > jr/ c ) ■ Br 2 (uu -> rr+7r 
4.0 x 10" 3 , 



where e = 0.96% is the efficiency determined from Monte Carlo and Br(uu — ► tt+tt tt°) = 0.89 [1]. 



VI. rj c -» </»</. 
A. Event Selection 



J/^i> — > 70(^ — » jK + K~K + K~ is used to study ry c — > </></>. The branching fraction of r\ c — > 00 in the PDG is 
mainly determined from the BES result Br{r\ c — > 00) = (2.5 ±0.5 ±0.9) x 10 -3 [13]. In order to do a systematic study 
of 7] c — ► FV, we also analyze the ry c — ► 00 decay, in order to remove the common systematic errors when calculating 
the ratios of reduced branching fractions BRn(r] c — > VV). 

The selection criteria for charged tracks and particle identification are the same as those in the r\ c — > p°p° 
analysis. At least three tracks must be identified as kaons. A 1C fit with a missing photon is performed to get better 
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momentum resolution. To select J/ip — > jcfxp candidates, \M K + K - — M^\ < 20 McV/c 2 , is required for both K + K~ 
invariant masses. 




2.8 2.85 2.9 2.95 3 

M 



KKKK 



FIG. 10: The r? c fit in the M u distribution. 



B. Results 

Fig. 10 shows the rj c fit to the distribution. The branching fraction of r] c <p(p is determined from 

N obs /e 



Br( Vc 



N Jfi> ■ Br(J/ij - 1Vc ) ■ Br^{4> -> K+K~) 
= (2.5 ±0.5) x 10~ 3 , 



where N obs — 76.1 is the fitted number of r\ c events in the distribution shown in Fig. 10, e = 15% is the 

efficiency determined from Monte Carlo simulation, and Br((j) — > K + K~) = 0.491 [1]. This result is consistent with 
the previously published BES result. 

VII. r\ c — > u:<p 

J/ip — > 77? c — > — > 77r + 7r~7T°i ; L + i ; r ~ is studied to obtain the rye — > w0 branching fraction. 



A. Event Selection 



The charged track selection criteria and particle identification are the same as in the r\ c — > p°p° analysis. At 
least two tracks must be identified as kaons. The photon selection criteria are different for r\ c — ► uxp to get higher 
efficiency because there are three photons in this channel. To reduce the number of spurious low energy photons 
produced by secondary hadronic interactions, photon candidates must have a minimum energy of 50 MeV and be 
outside a cone with a half-angle of 8° around any charged track. At least three good photons are required. 

To get improved momentum resolution and to remove backgrounds, events are kinematically fitted to 
the J/ip — ► 7777r + 7r - i ; r + if ~ (4C fit) and J/ip — ► 77r°7r + 7r~ K + K~ (5C fit) hypotheses, using all photons and 
7r + 7r~ K + K~ combinations. The combined probability is determined from particle identification probabilities of 
all tracks and the 4C fit probability to the J/ip — > jj-fTT + ir~ K + K~ hypothesis. Candidate events must satisfy 
xic(J/^ - U^+n-K+K-) < 25, X 2 C (J/^ ^n+n-K+K-) < 20, \M w+n - n o - M u \ < 0.03 GcV/c 2 , and 
\M K+K - - M | < 0.01 GeV/c 2 . 
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FIG. 11: The r\ c fit to the M 7T + 7I - 7T o K + K - distribution for candidate events. 



Results 



Figure 11 shows the result of fitting r\ c in the M w +^- v o K + K - distribution of J/ip — ► -fuj(f> candidate events. 
The shape of the r\ c used in the fit is obtained from i] c — > uj(p Monte Carlo simulation. The fitted number of r\ c events 
is 2.1 ± 6.1, which has large errors and includes decays other than J/ip — > "ftocp decay. So only an upper limit on 
Br(rj c — ► uj(f>) is determined. At the 90% confidence level, the upper limit on N Vc is N upper = 10. The upper limit on 
Br{f) c — > at the 90% C.L. is determined from 



jyupper j £ 



N 



Br(J/tP -> 7^) • Br{(j) -» K+K~) ■ Br{uj -> tt+tt-tt 



= 7.1 x 10~ 4 , 

where e = 4.2% is the efficiency determined from Monte Carlo. 



VIII. SYSTEMATIC ERRORS 



Many sources of systematic errors are considered. Systematic errors associated with efficiencies, such as from 
MDC tracking, particle identification, photon selection, and kinematic fitting, are determined by comparing ip(2S) 
and J/ip data with Monte Carlo simulation for very clean decay channels, such as ip(2S) — ► ir + n~J/ip. 

The MDC tracking efficiency has been measured using channels like J/ip — > AA and ip(2S) — > it + tt~ J/tp, 
J/ip — > n + fi + ■ It is found that the efficiency of the Monte Carlo simulation agrees with that of data within 1-2% per 
charged track. The total systematic error from the uncertainty of MDC tracking efficiency in our analysis is taken 
as 8%. The particle identification (PID) efficiency systematic error is calculated by comparing the efficiency of data 
with that of the Monte Carlo. According to the study of photon detection efficiency in J/ip — > pir [18], SIMBES 
simulates the photon detection efficiency in the full energy range within 1 to 3%. Here, 2% is taken as the systematic 
error in the detection efficiency for each photon. 

The systematic errors associated with the kinematic fit are caused by differences between data and simulated 
data in the momenta and the error matrices of charged tracks and the energies and the directions of neutral tracks. 
To check the consistency between data and Monte Carlo simulation, two channels, J/ip — > p°7r° and J/ip — > AA, are 
analyzed for 2-prong and 4-prong events, respectively. The systematic error for 4-prong events caused by the kinematic 
fit is determined to be 4%. For 4-prong events containing more than one photon, such as J/ip — > ■yq c — > -yuiui and 
J/ip — > 7^ c — ► juxp, 10% is estimated as the systematic error caused by the kinematic fit. 

The choice of different sideband regions can cause differences in Br(r) c — > VV). These differences are regarded 
as the systematic error associated with the choice of sidebands. Changing the fit range, binning, and background 
shape causes some differences in the result. The largest difference in Br(j] c — ► VV) caused by these changes is taken 
as the systematic error from these sources. 

The rj c mass resolutions used in the fits to determine the number of r\ c events are determined from Monte Carlo 
studies. Differences in r\ c mass resolutions between data and Monte Carlo also bring uncertainty to the determination 
of the branching fraction. The differences in Br(rj c — ► VV) caused by the uncertainty of r\ c mass resolution are taken 
as the systematic error caused by r\ c mass resolution. 

Many selection requirements are used to remove backgrounds. To estimate the systematic error caused by 
these, different selection requirements are used. The differences oiBr{r\ c — > VV) caused by the change of requirements 
are taken as systematic errors. 
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The number of J/tp events is (57.7 ± 2.7) x 10 6 determined from 4-prong events [10]. The uncertainty of the 
number of J/tp events is 4.7%. The branching fraction of J/tp — > -fr] c is (1.3 ± 0.4)%, according to the PDG [1]. This 
will contribute 30% to the systematic error of r/ c — > VV. 



TABLE III: Summary of systematic errors for r\ a — > VV (%). 



Sources 





T) c -> if* if* 


r\ c — ► cjij; 


r/ c — ► 00 


7/c — ► W0 


MDC tracking efficiency 


8 


8 


8 


8 


8 


Particle ID 


3 


3 


3 


3 


3 


Detective efficiency of photon 


2 


2 


10 




6 


kinematic fit 


4 


4 


10 


4 


10 


Different side-bands 


16 


12 








Fit range, binning and background shape 


3 


7 




11 




rjc mass resolution 


7 


8 




7 




Different Cuts 


18 


19 


14 


11 


32 


The number of J/tp events 


4.7 


4.7 


4.7 


4.7 


4.7 


Br(J/tp — > 7»7 C ) 


30 


30 


30 


30 


30 


Total 


40.7 


41 


37 


36 


46 



Table III lists all systematic error contributions and the total systematic errors. The Br(r) c — > VV) and 
Br(J/tp — > 7?7 C ) • Br(i] c — ► V^V") results including systematic errors are listed in Tables IV and V, respectively. Table 
V results do not have the large systematic error from Br (J/ip — > 77/c)- 

TABLE IV: Comparison of Br(rj c -> VV) values between BES and the PDG [1]. 



Br 


BES 


PDG 


Br(rj c pp) 
Sr(r7 c -► K'K*) 
Br(r) c — > ojo;) 
Br(r7 c — ► 00) 
Br(rj c — ► w0) 


(1.25 ±0.37 ±0.51) x 10 _;a 
(10.4 ±2.6 ±4.3) x 10~ 3 

< 6.3 x 10~ 3 

(2.5 ±0.5 ±0.9) x 10~ 3 

< 1.3 x 10" 3 


(2.6 ±0.9) x 1Q- Z 
(8.5 ±3.1) x 10~ 3 

< 3.1 x 10~ 3 
(2.6 ±0.9) x 10~ 3 



TABLE V: Summary of Br(J/tP 7 r7 c ) • Br(r) c -> VV). 



Br 


BES 


Br(J/tp — ► 7?7 C ) • Br(?7 c — > pp) 
Br(J/V> 7»7c) • Br(7? c JT7T) 
Br(J/tp — > 7?7 C ) • Br(7? c — > ujoj) 
Br(J/ip —> -)f) c ) ■ Br(rj c -» 00) 
Br(J/tp — > 7»7 C ) • Br(?) c — > u0) 


(1.6 ±0.6 ±0.4) x 10" 4 
(1.4 ±0.3 ±0.5) x 10~ 4 

< 6.7 x 10~ 5 

(3.3 ± 0.6 ± 0.6) x 10~ 5 

< 1.4 x 10" 5 



IX. SUMMARY 

The Br(rj c — > V^F) results are compared with PDG values in the Table IV. Table VI compares the ratio of 
branching fractions with expectations from SU(3) symmetry. The ratios arc more consistent with SU(3) symmetry 
than ratios based on PDG branching fractions, as shown in Table II. Our Br(r\ c — ► pp) result is consistent with the 
model prediction for T(rj c — > p°p°) of Ref. [3]. The upper limit on Br(r\ c — > w0) is also given. 
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TABLE VI: Check of SU(3) symmetry with new BES results. The common systematic error has been removed. 





BR(10- a ) 


BR R (W-*) 


Ratio of BRr 


SU(3) 


Br(rj c 


-> PP) 


12.5 ±4.9 


6.3 ± 2.4 


3.1 ± 1.2 


3 


Br{ Vc 


-> K*K*) 


10.4 ± 3.2 


6.1 ± 1.9 


3.0 ± 1.0 


4 


Br( Vc 


— > LOLO) 


< 6.3 


< 3.2 


< 1.6 


1 


Br( Vc 




2.5 ±0.7 


2.0 ±0.6 


1±0.3 


1 
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